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Abstract 

The momentum-space bosonization method of a Nambu and Jona-Lasinio type 
model with vector and axial-vector mesons is applied to irir scattering. Unlike the 
case in earlier published papers, we obtain the irn scattering amplitude using the 
linear and nonlinear realizations of chiral symmetry and fully taking into account the 
momentum dependence of meson vertices. We show the full physical equivalence 
between these two approaches. The chiral expansion procedure in this model is 
discussed in detail. Chiral expansions of the quark mass, pion mass and constant f n 
are obtained. The low-energy tttt phase shifts are compared to the available data. 
We also study the scalar form factor of the pion. 
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1 Introduction 



Coloured quarks and gluons are the main objects in the fundamental theory of the strong 
interactions, quantum chromodynamics (QCD), which are not directly observed in ex- 
periment. At intermediate and low energies the laws governing interactions of these 
fundamental particles bring about complex rearrangements inside a quark-gluon system 
which eventually result in two-quark and three-quark colourless formations — hadrons. A 
second transition appears at a comparable scale, namely the spontaneous chiral symme- 
try breakdown, leading to the pseudo-Goldstone boson octet. These phenomena cannot 
be consistently described within perturbation theory in the strong coupling constant a s 
which is no longer small at the typical hadronic energy scale. Therefore, it is mandatory 
to develop nonperturbative methods to solve these problems. 

There are few nonperturbative methods employed in particle physics. Noteworthy are 
lattice field theory, quantization on the light cone, and QCD sum rules. Despite all their 
advantages, these approaches have not yet lead to a unified description of confinement and 
spontaneous chiral symmetry breakdown (xSB). Both phenomena are supposed to arise 
from colour gauge interactions. Yet, the complex process of hadron formation remains 
obscure. This stimulates numerous attempts to construct models aimed at describing 
individual major features of the process. 

One can find explanation of spontaneous symmetry breakdown in the Nambu-Jona- 
Lasinio model, chiral quark model, approaches with stochastic vacuum fields, dual QCD, 
etc. A large class of such models involve path integral methods for calculation of Green 
functions that approximate QCD in the infrared region. The guiding principle is iterative 
integration first with respect to rapidly varying fields and then with respect to slowly 
varying ones, an idea borrowed from the general theory of phase transitions. According 
to this idea, slow modes mainly determine the asymptotic behaviour of the system at 
small energies and momenta (infrared asymptotics) and require thorough study in terms 
of collective variables. The contribution from "fast" fields is less substantial and can be 
taken into account by perturbation theory. The Nambu-Jona-Lasinio model (N JL) ^ , on 
which we concentrate below, is one of these schemes. The boundary between "slow" and 
"fast" fields is largely conventional being dictated by physical considerations. The chiral 
SU(3)l ® SU(3)r symmetry of the QCD lagrangian is broken at energies of the order of 
A x 5_b ~ 1 GeV if small quark masses are ignored. This is what allows one to relate A. x $b 
to the boundary between "large" and "small" momenta when looking into infrared QCD 
asymptotics from the symmetry point of view. 

Chiral symmetry breakdown is a fundamental issue because it is related to the problem 
of elementary particle masses and chiral symmetry lends itself to a controlable expansion 
(in energy) at energies below A x sb- At these energies the QCD vacuum becomes highly 
non-perturbative due to the condensation of quark-antiquark pairs as a consequence of 
the spontaneous symmetry violation. It is believed that any serious model of the strong 
interactions has to account for these facts. This is the case for the celebrated NJL model 
which was originally formulated not in terms of quark but rather nucleonic degrees of 
freedom. In what follows we will be concerned with a modern extension based on the 
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quark description. 

Assuming that the NJL model in the form that involves all essential QCD symmetries 
may be a reasonable approximation of QCD in the intermediate energy range, we develop 
a new theoretical method. It allows arbitrary N-point functions describing strong and 
electroweak interactions of basic meson states to be calculated through the effective four- 
quark interaction of the NJL model ||. Now there are two approaches to the problem. 
One is based on the classical treatment of four-fermion Lagrangians (the Hartree-Fock 
approximation plus Bethe - Salpeter type equations for bound states) f|-0. The other 
involves explicit separation of collective degrees of freedom in the functional integral of the 
theory (bosonization) P|-|p.4j|. In this framework one usually derives the Lagrangian based 



on a derivative expansion instead of directly calculating the Green functions. For technical 
reasons, such a derivative expansion can only be performed to the first few orders. We 
have constructed such a formal scheme that at the same time incorporates bosonization and 
has all advantages of the pure fermionic approach. It allows to combine the advantages 
of both of the abovementioned methods. The results of our more general method can be 
directly compared with the ones based on standard bosonization by performing a chiral 
expansion based on the nonlinear Bethe - Salpeter equations. 

Originally formulated in terms of boson variables linearly transforming with respect to 
the chiral group action, the method can be extended to nonlinear realization of symmetry. 
This is one of the significant advantages of the approach as compared with Hartree and 
Fock's classical technique, which is linear is this sense. Now chiral symmetry consequences 
are often investigated via nonlinear chiral Lagrangians. They are a series in energy (E), 
which is equivalent to normal expansion with an increasing number of derivatives. As a 
rule, investigations are mostly confined to the 0(E A ) approximation. Recently a program 
like this was used for the extended NJL model |TJ]]. In the present paper we show how 
one can generalize these results without any reference to an expansion in E 2 . 

Apart from pure theoretical developments to generalize the results of |TJj], we carry 
out explicit calculations. Our primary concern is the inr scattering amplitude. We shall 
calculate it both in the linear and nonlinear approaches using the standard Lagrangian 
of the extended NJL model with scalar, pseudoscalar, vector, and axial- vector four-quark 
interactions. We show that there is no difference between these two approaches. The 
expressions for the total amplitude A(s, t, u) are the same in both cases. It has to be 
so because the same Lagrangian is used. We consider this point in detail to show the 
self-consistency of our method. 

The paper is organized as follows. In the next section we introduce the Lagrangian and 
notation as well as review the main steps of the momentum-space bosonization technique. 
We then consider the nonlinear realization of chiral symmetry. We shall show that in 
this case equations for masses of bound states coincide with analogous equations in the 
linear approach. The expression for the weak pion decay constant /„- does not change 
either. All this allows us to get common rules for constructing chiral expansions in the 
linear and nonlinear approaches. In section 3 we derive the tttt scattering amplitudes and 
the main low-energy characteristics of the process in question. Theoretical derivations 
of the previous section will be the starting point. At first we shall use the model with 
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linear realization of chiral symmetry. In this section our calculations are an extension of 
earlier investigations [2, 15]. The corresponding amplitude includes the exchange of scalar 
and vector mesons as well as quark box diagrams. We give explicit expressions for this 
amplitude in terms of a few standard integrals. We obtain the low-energy expansion for 
the scattering amplitude and recover Weinberg's fl6j formula (to lowest order in the en- 
ergy expansion). Examining the model arising from classification of collective excitations 
based on nonlinear representations of the chiral group, we arrive at the same expression 
for the ittt scattering amplitude. We demonstrate this point by direct calculations and 
comparing the different contributions. In this part of the paper we not only generalize 



the known result of fll4| , we also prove the equivalence of this approach with the other 
known results. Our method looks like a bridge between different NJL models. Finally we 
present our results, comparing different variants of the NJL model with experiment and 
chiral perturbation results. In section 4 we calculate the scalar form factor of the pion. 
We pay particular attention to its relation to the value of the constituent quark mass. We 
conclude with a summary of our results. 



2 The theoretical background 



In Ref. [^] the extended NJL model was investigated and it was shown how to calculate 
Fourier transforms of iV-point Green functions by a simple method which was proposed in 
|2[ . This calculation was done in the linear realization of chiral symmetry. In this section 
we are developing the method to apply it to the extended NJL model with nonlinear 
realization of chiral symmetry. To describe collective excitations we shall use induced 
representations of the chiral group. Calculations of quark loop diagrams will involve 
renormalization of collective variables, as is typical of our approach. As a result, we are 
able to calculate amplitudes of physical processes to any accuracy in E 2 in models with 
nonlinear realization of chiral symmetry. 

2.1 The linear approach 



In this section we summarize previous results (see also paper 0]). The reader who is fa- 
miliar with the momentum-space bosonization method can use this section as a collection 
of our notation and main formulae to be found in this paper. The main idea of this method 
consists in the construction of special bosonic variables to be used for the description of 
the observable mesonic states. As a result, it extends the usual treatment of bosonized 
NJL models, which was formulated in || and developed in [9,10,12]. The standard ap- 
proach is essentially linked to the derivative expansion of the effective meson Lagrangian. 
The momentum-space bosonization method does not use this approximation. It involves 
the bosonization procedure and has all the advantages of the pure fermionic approach 
(Hartree-Fock plus Bethe-Salpeter approximation). Instead of step-by-step construction 
of the effective Lagrangian describing the dynamics of collective excitations we develop a 
method that allows direct calculation of amplitudes for particular physical processes. The 
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amplitudes derived accumulate the entire information on the process under investigation, 
just as if we had a total effective meson Lagrangian of the bosonized NJL model (in the 
one-loop approximation). 

Consider the extended SU{2) L ® SU(2) R NJL Lagrangian with the local four-quark 
interactions 



G 

L(q) = q{i^d^ - m)q + -y [(qr a q) 2 + (qi-f 5 T a q) 2 



G 



V 



{<lYr a q) 2 + {qi^r a qf 



(1) 



where q = (u, d) are coloured (N c = 3) current quark fields with current mass fn = 
diag(m u , rrid), r a = (r ,Tj), r = /, Tj {i = 1,2,3) are the Pauli matrices of the flavour 
group SU(2)f. The constants of the four-quark interactions are Gs for the scalar and 
pseudoscalar cases, Gy for the vector and the axial-vector cases. The current mass term 
explicitly breaks the SU(2)l <g> SU(2)r chiral symmetry of the Lagrangian ([!]). In what 
follows, we shall only consider the isospin symmetrical case fn u = frid — ffi- With boson 
fields introduced in the standard way, the Lagrangian takes the form 



L(q,a,n,v,a) = q (ry^ - m + a + i^it + + 757^^) Q 



°l + if 2 



2G.c 



v 2 4. n 2 

a _|_ iia 1 iia 



2G 



Here 



a 



OnTn 



TTnTn 



V 



(2) 



The vacuum expectation value of the 



scalar field do turns out to be different from zero (< do >^ 0). To obtain the physical 
field o"o with < 5~ >= one performs a field shift leading to a new quark mass m to be 
identified with the mass of the constituent quarks 



m 



On 



(3) 



where m is determined from the gap equation (see (|T3| ) below). 

Let us integrate out the quark fields in the generating functional associated with the 
Lagrangian (|2|). Evaluating the resulting quark determinant by a loop expansion one 
obtains 



L(cr, ii, v, a) 



—iTi In 

2 



1 + («7^ ~m) L (a + i7 5 7f + 7% + 757%) 



a 



7T 2 



2G 



+ 



S 



2G X 



(4) 



The NJL model belongs to the set of nonrenormalizable theories. Hence, to define it 
completely as an effective model, a regularization scheme must be specified to deal with 
the quark-loop integrals in harmony with general symmetry requirements. As a result, 
an additional parameter A appears, which characterizes the scale of the quark- ant iquark 
forces responsible for the dynamic chiral symmetry breaking. From the meson mass 
spectrum it is known that A ~ 1 GeV. Here, we shall make use of the a modified Pauli- 



Villars [18] regularization, which preserves gauge invariance and chiral symmetry. In this 
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form it was used in |19||-[p0[|. The Pauli-Villars cut-off A is introduced in the following 

way, 

e -im?z _^ R ^ = e -im^z [l _ (1 + lZ A 2 ) e - izA ^ 



m 2 e- im z 



iR'(z) = m 2 R{z) - izA 4 e- tz{A2+m2 \ 



(5) 
(6) 



where only one Pauli-Villars regulator has been introduced. In this case the expressions 
for the basic loop integrals Ii coincide with those obtained by the usual covariant cut- 
off scheme. Let us give our definitions for these integrals. To simplify the formulae we 
introduce the following notation 



A(p) 



9 9 ' 



d A q 



(7) 



Then we have 



h = iN c I d A qA{q) 

Up 2 



(4vr) 2 

iN c I d 4 qA(q)A(q+p) 
N, 



A 2 - m 2 In 1 + 



A 2 



1 '°°dz-, . |Vd-. 2 ) 



ir 2 / d y I — R ( z ) e 

l07T z jo jo z 



— R( z )e^ zp2{1 - y2) 



(8) 

(9) 
(10) 

(11) 
(12) 



Consider the first terms of the logarithm expansion in (^). From the requirement for 
the terms linear in a to vanish we get a modified gap equation 



UP 2 ) = ^ / dy{l-y 2 ) / —R{ 
62tc z jo jo z 

h(Pi,P2) = ~iN c J d A qA(q)A(q + Pl )A(q + p 2 ), 
h(pi,P2,Ps) = -iN c f d 4 qA(q)A(q + px)A(q + p 2 )A(q + p 3 ) 



m — m = 8mGsIi- 
The terms quadratic in the boson fields lead to the amplitudes 



U pp (p 2 ) 
Tl ss (p 2 ) 
U vv (p 2 ) 
U AA (p 2 ) 

U PA (p 2 ) 
U AP (p 2 ) 



8h-G s 1 +p 2 g- 2 (p 2 )\ <p+<pj,, 
8I l -G- s 1 + {p 2 -Am 2 )g~ 2 {p 2 ) 
g» v Gy l + A{p^ -g^p 2 )9v 2 (P 
'f (Gy l +Am 2 g- 2 {p 2 )) 
+ 4(pV - g^p 2 )gy 2 (p 2 )} e* A (p)e A (p) 

2img~ 2 (p 2 )p >1 e* A (p)ip P , 

-2img~ 2 (p 2 )p»E A (p)ip + P . 



(P)^(P), 



(13) 

(14) 
(15) 
(16) 

(17) 

(18) 
(19) 
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Here e^(p),e^(p) are the polarization vectors of the vector and axial- vector fields. We 
have introduced the symbols <pp = l and <p$ = 1 to explicitly show the pseudoscalar and 
scalar field contents of the pertinent two-point functions. The functions g(p 2 ) and gv{p 2 ) 
are determined by the following integrals 



<rV) = 4/ 2 (p 2 
9v\p 2 ) 



Let us diagonalize the quadratic form (14) + (17) 
fields 



4{P) 
(p) 



e£(p) - iP(p 2 )p^p, 
< A (P) + i(3(p 2 )p^t- 



This determines the function (3{p 2 ) 

P(P 2 



8ml 2 (p 2 



Gy 1 



lQm 2 l2(p 2 ) 



(20) 
(21) 

(19) by redefining the axial 

(22) 
(23) 

(24) 



Consequently, one has no more mixing between pseudoscalar and axial-vector fields. The 
self-energy of the pseudoscalar field takes the form 



U pp (p 2 ) = [8/x - Gs 1 + p 2 g~ 2 (p 2 ) (l - 2m(3(p 2 



(25) 



Now we can construct special boson variables that will describe the observed mesons. 
These field functions 00 correspond to bound quark - antiquark states and are derived 
via the following transformations 



where 



9AP 2 ) 



n a (p) 
a a (p) 

v a (p) 
a a (p) 

g(p 2 



Z^g^ip^ip), 

^ 1/2 <KpV(p), 

\z-^ 2 g v {p 2 )v a {p\ 

l -Z- l l 2 g v {p 2 )a a {p), 



'I - 2ml3{p 2 / 
The new bosonic fields have the self-energies 



g(p 2 



lQm 2 G v I 2 (p 2 ) 



(26) 
(27) 

(28) 
(29) 

(30) 



SahZ. 



-1 



2 111 
P -"V,<t(P 

p^Pv - g^v 



p~ 



m l,a(P 2 )]} ■ 



(31) 



2 Here and in the following we will use the common symbol 4> for the all set of meson fields: ir a , cr a , v a , a a . 
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The p 2 -dependent masses are equal to 



mlip 2 ) 
mlip 2 ) 

m 2 v (p 2 ) 
m 2 a (p 2 ) 



(G s 1 -8I 1 )gl(p 2 )=m(mG s ) 
1 - 2mf3(p 2 ) ml(p 2 ) + Am 2 , 

9 2 v(P 2 ) _ 3 



4Gv 



m v (p ) + 6m 



8G v J 2 (p 2 ) 

Hp 2 ) 



Hp 2 



(32) 
(33) 

(34) 
(35) 



These equations coincide with the conditions for appearance of quark - antiquark bound 
state as deduced in the pure fermion approach from analysis of the Bethe - Salpeter 
equations. 

The constants are determined by the requirement that the inverse meson field 
propagators 11* satisfy the normalization conditions 



ir-'V) 



p 



m 2 +0((p 



-fly 



around the physical mass points p 2 
values 



P - m v,a + 0{(P ~ m v,a 

mi, respectively. The conditions 



1 + 



ml [1 - 2mf3{m 2 K )} dl 2 {p 2 



Z v 

z n 



1 + 



1 + 



1 + 



m„ 



h(ml) 
- Am 2 dl 2 (p 2 ) 



dp 2 



himl) dp 2 
ml d.J 2 (p 2 



■M m2 v) 

2 

m„ 



dp 2 
dJ 2 {p 2 



Ji\ml 



dp 2 



6m 2 dl 2 (p 2 



Mm 2 



dp 2 



(36) 
lead to the 

(37) 
(38) 
(39) 
(40) 



In the following, when omitting an argument of a running coupling constant or a running 
mass, we always assume that its value is taken on the mass-shell of the corresponding 
particle. The symbol of this particle will be used for that. For example, on the pion 
mass-shell ml(p 2 = ml) = ml, j3{ml) = j3 n and so on. 

Using the expressions (|3ll) , one can obtain the two-point meson Green functions A^(p). 
For example, in the scalar and vector field case the relations 

U: b (p 2 )AUp 2 ) = S ac , K,(p 2 )A 



He 



(p) 



give 



K b {p 2 



$ahZ n 



p 2 — m 2 (p 2 ) 



Kr(p) 



5 ab Z v p^p u - g^ u m 2 v (p 2 
m 2 {p 2 ) p 2 — ml(p 2 ) 



(41) 



(42) 



This picture corresponds to the calculations in the framework of the pure fermionic NJL 
model where the Bethe-Salpeter equation sums an infinite class of fermion bubble dia- 
grams. 
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2.2 The non- linear approach 



Phenomenological meson fields with appropriate transformation properties under a non- 



linear action of the chiral group can be introduced as follows |22|. Let G be a 

continuous symmetry group of the initial Lagrangian and H a maximum subgroup of 
group G which leaves the vacuum invariant. Then an arbitrary transformation of the 
group G can be represented as G = K(()H(rj), where £, rj are the parameters determining 
the parametrization of the G group space. Acting from the left on the G group element 
by an arbitrary transformation of the same group G(g)K(QH(rj) = K((')H(i] r ), one can 
find out how the parameters ( and r\ are transformed under transformations of the group. 
It is essential that in this case the transformation for ( does not involve the parameters 
rj: £' = C'(C) <?)• Each parameter ( l is associated with a local Goldstone field ir l (x) so that 
the local fields ^(x) obey the transformation rule 

n*'(x)=C'(n l (x),g). (43) 

The non-linear transformation of the group G on the matter fields is constructed in the 
following way, 

Q^Q' = h(Tr,g)Q, (44) 

R—>R' = h(-K, g)Rh j (7r, g), (45) 

where Q, are the quark fields, R, the vector, axial- vector or scalar H multiplets, and 
h{ir,g)<=H. 

Let us consider the Lagrangian @. In this case G = SU {2)l®SU (2) r . The quark field 
q(x) can be represented as q(x) = ^(x) + q R (x), where qi{x) = Piq(x), qn(x) = P^g(x). 
The projection operators Pl,r are Pr l = (l±7 5 )/2. The fields qL,n{.x) transform linearly 
under action of chiral subgroups SU{2) L)R : 

<1l{x) -»• 9L{x)q L {x), q R (x) -»• g R (x)q R (x). (46) 



Let us introduce the notation a a + i^TCa = M a . Then we can write a + 2757? = MP R + 
M^P L . Now let us represent the complex 2x2 matrix M = M a r a as a product of the 
unitary matrix £ and the Hermitian matrix S 



M = £S£. (47) 

The matrix £ is parametrized by Goldstone fields, and its transformation law corresponds 
to the nonlinear transformation ( |4"3"D 

£ -> g L (x)^(rt, g L>R ) = h{it, g L , R )ig R {x). (48) 



The map £(71") : G/H — > G is thus the local section of the principal if-bundle G —>■ G/H, 
where h(jr, g) G H is a compensating H transformation which brings us back to our 
canonical choice for coset representative in the new coset specified by ir'. The exponential 
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parametrization = exp[iir / (2F)] corresponds to the choice of a normal coordinate 
system in the coset space G/H. New quark variables 



Qr = Ql = £ f <?L, Q = Qr + Ql 



(49) 



are transformed by the nonlinear representation of the group G, eq.(44), and can be 
used to describe the constituent quark fields in the approach under consideration. Let us 
rewrite the Lagrangian (§) as 



AGs 4GV 
Here we use the following notations: 



Q [ij^, + S - -(£ + 75 A) + 7 M (W^ - 75 W<- 
1 TrS 2 + -^Tr (wj^Wj^ + W^W^A . 



Q 



A = - £m£, 



To describe the vector and axial-vector * mesons we use new variables 



(50) 

(51) 
(52) 

(53) 
(54) 



We take the matrix S in the form S = m — m + s(x). Then the condition for the tadpole 
not to appear (after integration over quark variables) in the case of the scalar field s(x) 
will be the familiar gap equation ([13]). 

Among the terms quadratic in meson fields only the amplitudes 



H pp (p 2 ) = — 2m(m — m)Ii + (m — mj 2 p 2 I 2 (p 2 ) 



F 2 



U PA (p 2 ) 
U AP (p 2 ) 



m — m)p' l e : '*(p)l2(p 2 )(p l 



F 

8im , 



m — m)p^e l {p)l2{p 2 )^ 



(55) 
(56) 
(57) 



will be different than in the linear case (see (|T4|)-(19)). Let us consider a standard re- 
placement in a case like this 



(5S 



The condition for cancellation of nondiagonal terms (absence of the pseudoscalar - axial- 
vector transition) fixes the form of the function j3{p 2 ): 



hp 2 



8m(m — m)/ 2 (p 2 ) 
F[Gy l + 16m 2 I 2 {p 2 )) ' 
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(59) 



The self-energy of the pseudoscalar mode takes the form 



n Pi V 



4(to — to) 

Y 2 



2to 



TO — TO 



-h+V 2 h{p 2 ) 1 



2m(3(p 2 )F" 



TO — TO 



(60) 



Thus, in the nonlinear case only the pseudoscalar field 7r(p) will have a renormalization 
other than in the linear approach. The physical field 7r ph (p) should be introduced by the 
following replacement 

7r(p) = Z-V*&(pV*(p). (61) 

For other fields the transformations remain unchanged (see (27)-(29)). The function g n (p 2 ) 
is determined from (BUD 



F 2 [1 + 16to 2 GV/ 2 (p 2 )] = / Fg^ 
4(m — m) 2 l2(p 2 ) Vto — to 



(62) 



Hence it follows in particular that all equations for the masses of collective modes coincide 
with the ones in the linear case (|3^)-(35). This also applies to the equation of the pion 
mass, as seen from the chain of transformations 



ml(p 2 



to - m)gl(p 2 )Ii 



to[1 + 16m 2 G v I 2 (p 2 )} 
AmG s h{p 2 ) 



(63) 



where we employed the gap equation (|13|). The form factor f^ip 2 ) appearing in the vertex 
of the weak pion decay 7r — > £u# is easily found to be 



Up 2 



m 



(64) 



F[l + lQm 2 G v I 2 (p 2 )] VZ-gAP 2 ) 

The latter equality reveals that it coincides with a similar expression derived in the linear 
model [17|. Since the Bethe-Salpeter equation for the masses of collective states and the 
expression for the constant f n coincide in the two approaches under consideration, there 
exists a unified approach to construct chiral expansions. 



2.3 Chiral expansion 

At low energies the behaviour of scattering amplitudes or matrix elements for currents 
can be described in terms of Taylor expansions in powers of momenta. Yet, singularities, 
arising from the presence of light pseudoscalar particles in the theory, restrict the applica- 
bility of Taylor expansions. One has to take into account all these singularities to extend 
the divergence region for the series in momenta. This is possible because it is quite clear 
why the pion mass is small. The pion is a Goldstone boson and its mass is expressed 
in terms of current quarks which make it different from zero. Current quark masses are 
small and can be taken into account through perturbation theory. New combined Taylor 
expansion in momenta and current quark masses arising in this case form the basis of 
chiral perturbation theory [^3[ . 
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Let m and / be the values of the constituent quark mass m and the pion decay constant 
f n in the chiral limit where m = 0. In this case the pion mass is zero. The chiral series 
for this quantity begins with a term linear in m: 

The weak pion decay constant f n is (see eq.(64)) 

U = (66) 

Using (|30|) and (^Of), we arrive at 

fl = AZ- l 5m 2 I 2 {ml). (67) 

Hence, at fa — > we get 

o o o 2 o 

/ 2 = 45m J 2 , (68) 
which is non- vanishing in the chiral limit. Here and below we use the following notation 

h = h{0), 1 2= ]im h- (69) 

m— >0 

5 = 1- 2m/?^, 5= lim 5. (70) 
All these quantities are convenient to simplify the formulae. Note that the equality 

l-AG v f 2 =5 (71) 

is valid. We shall employ the gap equation (|T3|) to expand the constituent quark mass in 
a series in powers of the current quark mass,Q 

oo 

m = Y / c ( r ) m\ (72) 

Here 

c ( m ) =m= lim m, cf 1 ' = lim -ttt^t. (73) 

The mass m is a solution of the gap equation at m = 0. To find other coefficients of 
the series we differentiate equation flilf ) with respect to the current quark mass m. As a 
result, we get the following expression 

1 Til 

AG s M 2 (m) = . (74) 

m m 



3 We remark that only analytic terms in the current quark masses can appear within the Hartree 
approximation employed here. 
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We introduced the abbreviation M 2 (m) = 4m 2 J 2 . In the case of exact chiral symmetry 
with m = one can derive the value of M 2 (m) from equation (68). Bearing in mind these 
remarks, we get the following relation from the above equation: 

= i- (75) 

Other coefficients of the series are calculated through successive differentiation of the gap 
equation. For example, differentiating equation (74) we get 

, , , w mm' 1 m" , , 

8G S MM =— — -. 76 

m z m [m'Y 

On the other hand, relying on the definition of M one has 

MM' M 2 3h ( A 2 V 

mm m z 47r z yA^ + m^y 

In the chiral limit it follows from these two equations that 

(mh 2 / 2o ° 



(m) _ 3(4 ; ) 2 / m fei^ 
2m \ 2vr 2 / 



4 mJ = ^^|^-l). (78) 



For simplicity we use the symbol hi to denote the chiral limit 



m— >0 

Thus we obtain 



A 2 \ 

hl = \imhi= (^^2 • (79) 



sA 2 + m 



m =m 



o 2 o Q 4o2/ 2oo> 

4 m 32 m V Z7r / 



(80) 



Another important example is the chiral expansion for the pion mass. While in the 
case of the quark mass we employed the gap equation, here we need a pion mass equation 



([32]), which can be conveniently represented as 



m n -4mm— j / 2 ( m J = _. (81) 
It follows from the equation that this kind of expansion begins with a term proportional 

^ o 2 

to m. We have already called it m n (see (65)). Let us find the first few coefficients of the 
series 

oo 

ml = Y,<t^- (82) 
i=i 
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Obviously, Ci =m(Gsf 2 ) 1 . To find other coefficients of the series we represent equa- 
tion (81) as 

Gy 



.i=Q 



G, 



h(m 2 : 



1. 



3) 



Using the expansion (80) we isolate combinations to the same powers of m on the left- 
hand side of the equation and set them equal to zero. Thus we can calculate coefficients 
of the expansion (82). For example 



(~m n ) 
c 2 



(2f 2 G s y 



o 2o 

Tfl hi ° ° ° 

" L S (3 5-1) + (1- 2 8) 



2tt 2 / 2 



J4) 



S m *) 



o2 

8 



Am{2PG s ) 



S5-3 + 



o 2o 

m hi 
2tt 2 / 2 



o o o2 2 o 

5(2-105+215 )m hi 



2vr 2 / 2 



+75-205 



2 4 4(155 2 -105+2)m 2 
» 5(A 2 + m) 



J5) 



Below we give the results of similar calculations for the pion decay constant f n on the 
basis of the expression 



In this case the corresponding coefficients of the chiral series are 



Jf*) 



(86) 



If) 



8 



AmfG s 



o 2 o o ' 

3 m hi8 



c 2 



o3 o 

35 (2-5) 



o2 o 

5 hi 



2f (AmfGs) 2 2/(4tt/ 2 G s ) 



2 

35 



A 2 + 2m 
A 2 +m 



o 2 o o 

9ra fa5 



(2-5) 



o 2 o o 

m hi8 



o 2 

A 2 + 3m 



(2tt/) 2 5 (A 2 + ^ 2 ) 



59) 



Concluding the section we point out that chiral expansions derived here for the main 
meson characteristics prove to be helpful in establishing correspondence between the re- 
sults obtained here and the known low-energy theorems of current algebra. We shall use 
them to analyze -an scattering and to consider the scalar radius of the pion. 
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3 7T7T- scattering 



The formal scheme described in the previous section gives the possibility of evaluating 
any mesonic N-point function through the parameters of the model: A,m, G$, Gy. Let us 
use it for the evaluation of the 7T7r-scattering amplitude. We have done that in the same 
approach earlier ||, but without considering the spin-one mesons. Now we would like to 
explore the role of vector and axial-vector particles in this process. 

The model with linearly realized chiral symmetry was already employed for this purpose 
||15|| . However, the results of [15 should be considered only as a first approximation to 



this problem. The it — a\ mixing was neglected there. This breaks chiral symmetry, 
which is recovered only if the constituent quark mass m goes to infinity. In our case we 
exactly reproduce Weinberg's result at the level E 2 and we also evaluate all higher order 
corrections in E 2 . 

The extended NJL model with nonlinearly realized chiral symmetry has not been used 
to study low-energy tttt scattering so far, though some general conclusions can be drawn 



from the results of |]I4| . Chiral expansion parameters Li were calculated in the model, 
the standard heat kernel method being used to get the first terms of the effective meson 
Lagrangian (up-to-and-including the 0(E 4 ) terms). We shall advance farther and find 
a general (to any accuracy in E 2 ) form of the tctt scattering amplitude in the principal 
approximation in 1/N C . 

In our calculations we use the conventional Mandelstam variables: 

s = {qi + g 2 ) 2 , t=( qi - q 3 ) 2 , u = {q x - q 4 ) 2 (90) 

for the scattering process vr a (gi) + 7r fe (g 2 ) — » ?i" c (<?3) +vr d (g 4 ). The nn scattering amplitude 
T a b;cd has the following isotopic structure 

T a b;cd( S > *! u ) = 8ab8 c dA(s, t, u) + 5 ac 5 bd A(t, s, u) + 5ad5 cb A(u, t, s). (91) 

It follows that the amplitudes with definite isospin are 

T°(s, t, u) = A(t, s, u) + A(u, t, s) + 3A(s, t, u), 
T 1 (s, t, u) = A(t, s, u) — A(u, t, s), 

T 2 {s, t, u) = A{t, s, u) + A(u, t, s). (92) 

Meson tree diagrams corresponding to the amplitude A(s, t, u) are of the same form both 
in the linear and nonlinear approach (see Fig.l). Only the internal structure of quark- 
loop-based meson vertices will be different. 

3.1 The amplitude (linear approach) 

Let us calculate the amplitude A(s,t,u) in the linear approach. The vector meson and 
scalar meson exchange diagrams of Fig.l include the triangular vertices p — » tttt and 
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a — > 7T7T. We calculate these vertices in accordance with the diagrams of Fig.2, taking 
into account the nai mixing. The p®(p) — > 7r 6 (p 1 )7r c (p 2 ) amplitude is equal to 



M, 



/miT ^ 



Tr (r a [n , t c ] - ) (pi - Vi ) % (p) /p™ (pi , P2) , 



where 



(93) 
(94) 



The function F(pi,p 2 ) has the form (for on-shell pions) 



^) = t1- 



+ 



m 2 {p 2 ) ' p 2 — 4m 2 



A hK-PuPi) 
h{ml) 



(P 2 - 2ml) 



' hip 2 ) 

MO 



- 1 



+2m 



(95) 



Here p = P1+P2, with pi,p 2 being the pion momenta. This function obeys the requirement 
of universality of electromagnetic interactions, F(0) = 1. To see this one has to use the 
equality 

,dh(p 2 ) 



J 2 (0) - J 2 (m 2 ) - mih(- Pl ,p 2 )\ p2=0 = 2m 2 - 
The a a (p) — > 7r b (pi)n c (p 2 ) amplitude has the form 



dp 2 



(96) 



M an7T (p 1 ,p 2 ) = -Tr (r a [T b ,r c } + ) f an7r (p 1 ,p 2 ). 



(97) 



The vertex function fwitiPiiPi) on the pion mass shell is 



farnv ip ) 



16mglg(p 2 



p z 



Am 2 



[1-5 2 ) 



m 2 5 



W) + ^(i-5)/ 2 «) 



+ -(p 2 -2ml)I 3 (- Pl ,p 2 ) . 



(98) 



Using these expressions we obtain the contributions of the meson-exchange diagrams 
in Fig.l to the amplitude A(s,t,u), 



A p (s,t,u) = 4GV 



(s-u)F 2 (t) (s-t)F 2 (u) 



l-\G v tJ 2 {t) \-\G v uJ 2 {u) 

Zcrfa7Tn{ S ) 



A„(s,t,u) = 



(99) 



(100) 



m 2 (s) — s 

Now let us consider the remaining diagrams, the boxes. Their total number is 16. 
In Fig. 3 we display only structurally different diagrams. Their calculation is quite 
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cumbersome. Here we collect only the final result, 

s 



2^2 

m 2 



A rest (s,t,u) = AgtZ- 2 ^A[sl3%{ml) - I 2 {s)} + (1-5 

x [I 2 (s) - I 2 {ml)\ + lft[(s - u)tJ 2 {t) + (s- t)uJ 2 (u)} 

- 6 2 [4(s - 2ml)h( qi , -q 2 ) - ^C(s, t, u) + 5 2 B(s, t, u)}} , (101) 



where 



B(s, t, u) = (2m 4 . - us)I A (q 1 , q x + q 2 , q 4 ) + (2m 4 - ts)I A (q 1 , q x + q 2 , q 3 ) 

-(2m*-ut)h(qi, qi -q 3 ,q 4 )). (102) 

C i s ^) = " 2m ')[ J 2W " J 2K)] + 2^J 3 ( gi , g 3 )} + (*<-> u). (103) 

Let us consider the low-energy expansion of the total 7T7r-scattering amplitude 

A(s, t, u) = A p (s } t, u) + A a (s } t, u) + A rest (s, t, u). (104) 



The self-consistent procedure to do that is chiral perturbation theory p3[ . The chiral 
expansion in powers of external momenta and current quark masses in the case of the 
scattering amplitude has to lead to Weinberg's celebrated theorem [pi at the E 2 level, 

o 2 
e — 777 

A(s,t,u) = —^ + 0(E 4 ), (105) 

o 2 o 4 

where 0{E i ) is the short form for the terms of order q A , rn n q 2 , and higher. In the 
case at hand the simplest way to get this result is to consider first of all the momentum 
expansion and only partly use the quark mass expansion. We shall do a chiral expansion 
for the pion mass, m^, and the constituent quark mass, m, only at the last stage. In 
accordance with these remarks one obtains 

A >="4 ( 3 -^)Hl-S)I^O { E% (106) 



.4 



rest 



A a = ^$Ul 2 + ^- s [25 2 -l) ,niH U) - 3) 



Z 2 



87r 2 m 2 



AI 2 (l-(3 2 s) + 



s - 35 2 (s - 2ml 
/,J .s - 3^(.s - 2///: 



0{E A ). 



]ir 2 m 2 



+ 0(E 4 



(107) 



(10£ 



Some remarks are useful here. We use momentum expansions of integrals his) = I 2 + 
(shi) I (327r 2 m 2 ) + . . . , I-i{qi,q2) = — 3/ii/(327r 2 m 2 ) + . . .. One can see that not only 
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divergent terms (~ Iq) contribute at this level. There are finite contributions which are 
proportional to h%. Full cancellation of these terms in the full amplitude 



A(s, t, u) 



zi 



s — m„ 



+ 0(E 4 



(109) 



is a good test of self-consistency of our approach. In order to preserve chiral symmetry 
they must start only at the next-to-leading order. Noting that 



77T 



1 

P 



0{E 2 



110) 



we have full agreement of our amplitude with the low-energy theorem (105). 

Now let us consider the E A correction to Weinberg's result. After some lengthy calcu- 
lations one obtains 



Here 



12vr 2 / 2 ° 2 



m 



96vr 2 / 4 



o 2o 

m hi 



2Ii(s-2m I ) 2 + h s 2 + (t 



u 



+ 



o3 4 o 2 

95 m hi 



7T 2 / 2 (2*/)' 



o4 

35 h2 



o 2o 

m § 



and 



hi- 



o2 

:i-s) 



o4 

95 h 2 



o2 o2 

65 (1-5 ) 



47T 2 / 2 

o 2o 

m 5 



o2 



A 2 + 3 m 

^2~ 

A 2 + m 



112) 



113) 



(114) 



Let us compare our result (111) with the known estimations of the 0(E 4 ) terms. One 
can see that the structure of the term A^(s, t, u) fully agree with the early result of [2~1|. 

_ _ o 

There parameters l\ and h, in the limiting case, when Gy — > 0, i.e. 5^ 1, go over into the 



known result PD]| , 0. In the case of 577(3) symmetry this set of parameters are known 
as Lj. They have been obtained in the paper [[14]] on the basis of S77(3) ® SU(3) NJL 
Lagrangian with spin-one mesons. The constants Li specify the general effective meson 
Lagrangian at order of E A . One should, however, not compare them directly to those 
analyzed by Gasser and Leutwyler since their analysis includes the effect of meson loopsfj] 



If one considers the SU (2) limit of the effective meson Lagrangian from [14[] and compare 
it with the known SU{2) effective meson Lagrangian from p0[ one can relate the S77(3) 
low-energy parameters Li with the 577(2) constants k- In particular we have 



L 2 = 2Li 



h l\ — I2 

192vr 2 ' 3 ~ 1927T 2 ' 

And thus leads to a scale-dependence not present in the NJL model. 



115) 
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Or 



L2 — 2L\ 



64vr 2 



27T 2 / 2 

o o 2 
35m 



5 foH 



o2 o2 Q 2 o 

'1-5 f + 25 (1-5 ) /11 



116) 



64tt 5 



o4o 27T 2 f 2 o2 o 2 2 o 

"45 /12— S ) 2 + 25 (45 -3)/u 
5m 



o3 / O n 2 O V ^ 

47r 2 / 2 5 / x 35m 



o 2 

m 



4tt 2 / 2 
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To compare our formulae with similar expressions derived in [fL4|| , we point out the fol- 
lowing rules of correspondence: 



r(o,x) 



4tt 2 / 2 

o o 2 
35m 



, r(i,x)~/u, r(2,x)~/?,2, 9a ^5 ■ 
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Here, we give the notation of fl4| on the left side and our equivalents on the right side. 
Now it is easy to make sure that this part of our calculations yields results which agree 
fully with earlier estimates. 

Let us calculate scattering lengths aj and effective range parameters b\ up to and 
including 0(E 4 ) terms. The result is given in the Appendix since it is not new. The 
subsequent numerical estimates will reveal that the E 4 approximation of the extended 
NJL model is a very good approximation to the Hartree solution. 



3.2 The amplitude (non-linear approach) 

Let us calculate the tttt scattering amplitude in the nonlinear approach. Though we used 
the same initial Lagrangian, the scheme developed differs from the approach with linearly 
realized chiral symmetry. The difference arises from variation in the structure of vertices 
describing the interaction of Goldstone particles with quarks. Note that the structure of 
similar vertices for scalar, vector, and axial- vector fields remains unchanged. As before, we 
begin with the amplitude p^(p) — > ^ b (pi)TT c (p2) (see (93)). To distinguish linear results 
from their nonlinear counterparts, we shall use symbols with a tilde above them. For 



example, after calculation of the diagrams in Fig. 
(95)) we get a new expression F{p 2 ): 



4, instead of the function F(p ) (see 



hp 2 

(mg A - m) 
p 2 — Ami 



4fi 



Z w (m — m) 



P 
12 



(1 - QaJMp 2 ) ~ m(mg A - m)I 2 {m\ 



(p 2 - 2m 2 ) (l 2 (p 2 ) - J 2 (m 2 )) + 2m%(-p uP2 ) 



(119) 



The constant g A is equal to 1 — 2/3 n F. It can be related to the known constant 5: (mg A — 
m) = (m — m)5. If m, p 2 — > 0, then F(p 2 ) — > 0. This off-shell behaviour of the form 
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factor essentially differs from the behaviour of F(p 2 ), where F(0) = 1. The contribution 
to the 7T7T scattering amplitude from the p meson exchange diagram will be described by 
formula (99), where one should make a change F(p 2 ) —>■ F(p 2 ). We call this contribution 
A p (s,t,u). For the given amplitude, owing to the above property of the function F(p 2 ), 
the chiral expansion begins only with the terms 0(E 6 ). 

The vertex function for the two-pion decay of a scalar particle^ will also differ from the 
linear one. It is described by the second group of diagrams in Fig. 4. 



fa-Kir (p 



4g(p 2 )gl 



+ 



7 -y-, -W {C 1 ~ 2g A )(mg A - m)m 2 J 2 (m 2 n ) 

mg 2 A + m{l — 2g A ) p 2 I 2 (p 2 ) — 4mrfi(m — m)I 2 (p 2 ) 



+2m(mg A - m) 2 (p 2 - 2mJ)/ 3 (-p 1 ,p 2 )} . 



(120) 



The form factor f a7Tn (p 2 ) also tends to zero at m, p 2 — > 0. Changing f an7T for f a7Tn in (100) 
we get the contribution to the tttt scattering amplitude from the exchange of the given 
scalar. We denote it by A a (s,t, u). The chiral expansion of A a (s,t, u) begins with terms 
of order 0{E A ). 

Now let us consider the remaining group of diagrams displayed in Fig. 5. Here we do 
not show diagrams resulting from elimination of ita\ mixing. These diagrams are easily 
derived from the given ones by the corresponding inserts in vertices £ with one pion at the 
end. Without going into detail, we shall give the total contribution to the tot scattering 
amplitude from all the diagrams taken together, 



A 



rest 



S,t,U) 



Z 2 J 



m 



m) 



{{mg A - 



m) 



im 



3m) (s — m 2 ^] 



+2m(l - 2g A )ml / 2 (mJ) + 



2 \2 



;i - g\ ) 
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[{s - u)tJ 2 {t) + {s - t)uJ 2 {u)} 



+m 



m(l — g A ) 2 s + 2g A {mg A — m)s — 4m(m — m) 2 his) 



+g\{mg A - m)(mg A - 3m) s I 2 (s) - h{m\ 
+ {mg A — m) 2 4m(m — m)(s — 2m 2 )I 3 (qi, —q 2 ) 
-(! - g 2 A) c ( s i *> u ) ~ ( m gA - m) 2 B(s, t, u)] } . 
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Here we use the notation (102) and (103) adopted earlier. 

The expressions for the separate amplitudes A p ^ rest (s, t, u) in the linear and nonlinear 
approaches are different. Before we will consider the total amplitude A(s,t,u) = A p + 
A a + A rest let us show that its chiral expansion coincide up to and including the terms 
0(E 4 ) with the corresponding result in the linear approach. Indeed, in the main 0(E 2 ) 



5 Similarly to the linear case, we shall use the symbol a for this particle, though in the Lagrangian the 
symbol s was used. 
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approximation we get Weinberg's result which follows from the term 



Z 2 (m — m) 4 



{rag a — m)(m — 3m) (s — m^^im^) 



o 2 



p 



(122) 



In the 0(E A ) approximation the contribution comes from only two amplitudes, A a (s, t, u) 
and A rest (s, t, u). As was established, p meson exchange diagrams lead to an amplitude 
beginning with the term of the order 0(E & ). Here we get 



A 



(4) 
rest 



>mm 



f 4 



° ^ o2 o o o 2 

2mm — s 5 -5 (1 -2S)m 7r 



Am f 2 



Smm 



lis -2ml) I 1 



o 2 o o 

3 m Sh i 
4tt 2 / 2 



°h+^-k(s-2ml- 
2 16tt 2 V 

i 2 

o 2 ' 



o o2 o o 2 

2mm - sS -8 {l — 2S)m 7r 



o 



(123) 



(1-5 ) 2 



o 2 o 



6/ 4 



Am p 



2s°S +(1-45) ml 



o2 o 

5 hi 



All I , o 2 o2 2 



3s(l-(J)+2m w (35-2) 



4 
5 /i2 

"8vr 2 / 4 



+(!-? - + (s - *)(u+m*)] 



2m n +ut — s(t + it) 



(124) 



Summing the expressions, we arrive at a result coinciding excately with the result of simi- 
lar linear calculations (111). Noteworthy is that contributions proportional to m are fully 
cancelled by the summation. Thus, despite all the differences in the two approaches, they 
are equivalent up-to-and-including the terms of the order 0(E 4 ) in the chiral expansion. 



3.3 Equivalence of the linear and the non-linear approach 



Haag's theorem |25| in axiomatic field theory states the independence of the S-matrix 
elements on mass-shell from the choice of interpolating fields. In the framework of the 
Lagrangian approach the same result exists [P6|j . m our case it means in particular that 
if we did all correctly the total amplitude A(s, t, u) should be the same in both cases. To 
demonstrate it let us first consider the pun form factors E(p 2 ) and E(p 2 ) (see (95) and 
(119)). One then obtains 

F{p 2 ) - F{p 2 ) = ±- 



m 



m — m 



V 

m 2 (p 2 



- 1 



(125) 
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It is clear now that on the p mass shell these form factors coincide with each other. We 
have a similar behaviour for f an7V and / CT7r7r (see expressions (98) and (120)), 



fairirij) ) faTnr(p ) 



Z^\fZl{m — m) 



p - m ff (p ) 



Again at = m 2 one gets zero for this difference. 
We are ready now to compare the amplitudes 

A Q (s, t, u) = A a (s, t, u) - A a (s, t, u), 

where a = p, a, rest. The result can be written as 



A p (s,t,u) 



4gW - S) 

Z 2 m(m — m) 



2m — m 



m — m 



(3s - Aml)I 2 {ml) 



-D(s,t,u) + 25C(s,t,u)]. 



(126) 



(127) 



(128) 



A a (s,t,u) 



4^ 



Z 2 (m — m) 2 

s(l-5 2 



[ml5I 2 (ml) + (Am 2 - s)I 2 (s) 



—Am(m — m) 2 1 — 

I \ Am 2 

+5 2 (s-2m 2 7r )h(q 1 ,-q 2 )}}. 



m 



h(s) + -^5(l-5)I 2 (mi) 



m 



(129) 



A res t(s, t, u) 



s + 4m(m - 2m) 2s(l - 5 2 



m 2 5 



m(m — m) 2 
Am5 2 



(m — my 
3s 



m(m — m) 



Am(l —8) — m H -(1 — S)(m — 2m) 

m^ 



h(s) 
h(ml) 



+ 



[m — m 
m(m — m) 



s -2ml)I 3 (q 1 ,-q 2 ) 
[D(s,t,u) -25C(s,t,u)}\ . 



Here the notation 



65m 2 



m — m 



[(s-u)U 2 (t) + (s-t)uJ 2 (u)} 



has been used. Summing (128)- (130) it is easy to see that 

A p (s, t, u) + A a (s, t, u) + A rest (s, t, u) = 0. 
This is what we wanted to show. 



(130) 



(131) 



(132) 
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3.4 The full result 



Now we can go on and estimate the role of higher order contributions in the amplitude 
obtained. For direct comparison with the empirical numbers we fix parameters in order 
to obtain the pion decay constant, /„■, and the pion mass, m^, close to their physical 
values, f n ~ 93 MeV and m^ ~ 139 MeV, respectively. We take two sets of parameters, 
set / with a constituent quark mass m ~ m p /2, set II with a rather low quark mass. 
In the first case we get a bound state for the p - meson, in the other case the p - meson 
mass lies above the threshold for production of a pair of 'free' quarks, a problem related 
to the absence of confinement of the model. Nevertheless, as will be discussed in the 
next section, a better description of the scalar form factor of the pion is only achieved for 
rather low constituent quark masses, so that, with due care, it is of interest to study pion 
observables for this case too. The model parameter Gy is obtained by fixing the mass of 
the p meson, m p ~ 770 MeV, as long as m p < 2m. For the low constituent mass case, 
the p is embedded deeply in the qq continuum. To avoid the complications of defining 



this isovector- vector state under such circumstances [27], we choose to fit the scattering 



length a \ to fix Gy. This is an unambiguous and rather simple procedure. As a result, 
the four parameters of the model for set / have the values: G$ = 9.41 GeV -2 , Gy = 
11.29 GeV -2 , m = 390 MeV (m = 3.9 MeV) and A = 1 GeV. With these parameters, 
we find f n = 92 MeV, m w = 139 MeV, m p = 770 MeV, 5 = 0.62. For set II we obtain 
G s = 1.083 GeV -2 , Gy = 8.8 GeV -2 , m = 200 MeV (m = 1.0 MeV), A = 2.5 GeV, 
/„- = 92.7 MeV, m n = 139 MeV, 5 = 0.69. For obvious reasons, phase shifts are always 
presented for -^s < 2m. We also calculate the set of parameters in the limiting case 
Gy — > or S — > 1, already determined in [fj. Here we use Gs = 7.74 GeV -2 , A = 1 GeV 
and m = 242 MeV (m = 5.5 MeV). Therefore one can always compare the predictions of 
the NJL and extended NJL models. Apart from that, we need values of the main physical 
quantities in the chiral limit m — > 0. In this case we obtain for set /, m= 382 MeV, / = 
90.9 MeV, 5= 0.628 and for set II, m= 183 MeV, / = 88.4 MeV, 5= 0.73 The leading 
term of the chiral expansion for the pion mass is m 7T = 138.42 MeV in / and m n = 141 
MeV in set II. 

First, we give the values for the constants L2 and L3: we have for set /, L2 = 1.2 • 10 -3 , 
L 3 = -3.2 • 10 -3 and for set II, L 2 = 2.0 • 10 -3 , L 3 = -2.2 ■ 10 -3 . For comparison, 
we give the scale-dependent empirical values L 2 {rn p ) = 1.2 • 10" 3 , L 3 (m p ) = —3.6 • 10 -3 . 
However such a comparison has to be taken cum grano salis since the L t calculated 
within the Hartree approximation are simple c-numbers. Whereas set / lead to values 
closer to the experimental ones, set // is compatible with the results of |TJJ]. Second, we 



consider the im threshold parameters. The results of our calculations are presented in 
the Table. The first and second columns list the numbers based on the chiral expansion 
(E A approximation), for sets / and II. The formulae used are given in the Appendix. 
The third and fourth columns show the results of the exact calculations on the basis of 
the full amplitudes A(s,t,u), for I and II. We compare them with the results from the 
standard NJL approach without vector mesons 0], within the current algebra |TH], and 



with the experimental data [28 
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Table: The tttt scattering lengths and effective ranges in the extended NJL model in comparison 
with the same calculations in the original NJL model, without spin-one mesons and the soft 
meson theorems (SMT) |L6|]. The experimental data are taken from [28|. The '*' denotes an 
input quantity. 



From this comparison one infers the following. All scattering lengths and range pa- 
rameters are mainly determined by the 0(E 4 ) approximation. The effect of higher-order 
terms is noticeable only in the D-waves. If we compare the results of numerical calcu- 
lations with our previous estimations of low-energy tttt scattering parameters within the 
NJL model without vector mesons, we notice that the agreement with the experimental 
data is of comparable quality (for set II) but poorer for set I. This confirms the obser- 
vation of Ref. |2{| that within the Hartree approximation one needs to work with a small 
constituent quark mass. Similar findings were obtained in Ref.|[L4||- Alternatively, one 
could think of going beyond the Hartree approximation and, in particular, include pion 
loops. This is a difficult problem which deserves further studies ||30|| . 

In Fig.6a,b we show the S and P-wave phase shifts 5q, 5q and S\ for set I in compar- 
ison with the available data (-^/i < 700 MeV). The calculated phases are in reasonable 
agreement with the data. For set II, we are confined to ^Js < 400 MeV and therefore 
only show the S-wave phases in Fig. 6c. New data in this region will come once DA$NE 
is operating and Ku decays have been analyzed. 



4 The scalar form factor of the pion 

Another quantity closely related to tttt scattering is the scalar form factor of the pion, 
denned by 

< ttV) I m(uu + dd) | n\p) >= S^T^t) (133) 
where t = (p — p') 2 is the square of the invariant four-momentum transfer. In the frame- 



work of the original NJL model this quantity was previously determined |34[] and agreed 
quite well with the empirical scalar form factor of the pion for low momentum transfers at 
a fairly small constituent quark mass of 241.8 MeV. In particular, the scalar pion radius 
was demonstrated to impose powerful constraints on the parameters of the NJL model. 
We anticipate that with the large constituent mass as in set I, the scalar radius will be 
largely underestimated. The relevant Feynman diagrams to calculate this quantity are 
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depicted in Fig. 7. There is a sum of two contributions to the scalar pion form factor in the 
NJL model: one is the direct coupling of the operator m(uu + dd) (double-dashed line) to 
the two pions via a quark triangle (we call it bare coupling), and the other corresponds 
to rescattering of quarks into the scalar meson a, which in turn couples to the pions. 
In Fig. 7a we show these contributions, with the pion legs amputated. After solving the 
Bethe-Salpeter equation for this vertex, T, one obtains the scalar pion form factor by 
attaching the pion legs, Fig. 7b. With eqs. (20) and (98), the result for the scalar form 
factor of the pion T n (t) can be cast in the form 



where 



JC(t) 



I\r(t) 



1 + 



J s {t)g\t) 
ml(t) - t 



771 

hit) + mlPJI 2 (ml) + -5\t - 2ml)I 3 (p,pi) 



(134) 



(135) 



and J7s(t) = 8ii + 4(£ — 4m 2 ) hit) is the fundamental quark bubble in the scalar channel. 
For the same parameter sets used in the evaluation of the scattering parameters, we 
obtain r 7r (0) = 1.007m 2 , (set I), and T w (0) = 0.946m 2 , (set II), consistent with the xPT 
prediction |23| . By expanding in powers of t we extract the scalar mean square radius of 
the pion 

Wi =i+ l <r2Kt+o{t2) - (136) 

Its numerical value < r 2 = 0.043 fm 2 calculated with the parameter set / is one order of 



magnitude smaller than the empirical value < r >^= (0.55±0.15)fm [29 extracted from 



phase shift analyses |3~5| . For set II we increase the mean scalar radius < r 2 >*= 0.53 fm 
in good agreement with the empirical number. To understand these numbers in more 
detail, we investigate the chiral expansion of the scalar form factor up-to-and-including 
terms of order E A , 



r w (o) 
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(137) 



In the leading order of the chiral expansion we get exactly the current algebra result 

o 2 



r 7r (0) =m 7r ||36|| , as expected, since all symmetries are respected in the evaluation of the 

o 2 

form factor. The E contribution leads to r 7r (0) = 1.005 m^ for set / and r 7r (0) 

o 2 

0.943 m n in set //. For the scalar radius of the pion we have 



o2 
o 2 

2m 



o 2 o o 

3m Shi 



(138) 



which is ~ 0.057 fm 2 in set I and ~ 0.70 fm z in set II, i.e. in this quantity one observes 
large higher order effects. This is not unexpected from previous calculations in chiral 
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perturbation theory beyond one loop []35 |. These higher order effects are essentially ac- 



counted for by the low constituent quark mass in agreement with the findings of Ref. |2P 
We stress again that alternatively one might want to consider pion loop effects, which, 
however, goes beyond the scope of the present manuscript. 



5 Summary and conclusions 

In the present paper we have developped a method which allows for the calculation of Ap- 
point functions within the extended Nambu-Jona-Lasinio model with nonlinear realization 



of chiral symmetry [|14]]. This extends the ideas described in 0. Using the path integral 
technique and the Hartree approximation, we work in momentum space to circumvent the 
standard way of employing the heat kernel expansion for constructing an effective meson 
Lagrangian. With the heat kernel method, one manages to get only the first few terms 
of the Lagrangian (in general the 0(E 4 ) approximation) and faces growing difficulties 
when calculating terms of yet higher orders. Examining the one-loop approximation for 
the effective action in momentum space, we find transformations of collective variables 
which permit to calculate these higher order terms in a straightforward manner. The two- 
point functions describing propagation of collective excitations lead to results coinciding 
with those obtained from analysis of Bethe-Salpeter equations on bound quark- ant iquark 
states for two-particle Green functions of quark fields. 

Using equations for the mass of constituent quarks (gap equation) and the pion mass, 
we construct chiral expansions for these quantities and the constant f n . We have also 
calculated the tttt scattering amplitude. This program is applied both to the linear and 
nonlinear description of transformation properties of physical fields. We have demon- 
strated explicitely that the tttt scattering amplitude A(s, t, u) obtained in both approaches 
is completely equivalent. We have carried out the chiral expansion for this amplitude and 
calculated the pertinent scattering lengths and range parameters. In addition we have 
obtained the corresponding low-energy tctt phase shifts. All these calculations have been 
performed for two sets of parameters, the first and the second one having a large (~ 400 
MeV) and a small (~ 200 MeV) constituent mass, respectively. In general, the conclusion 
is that the corrections derived by considering terms of the effective meson Lagrangian 
with higher derivatives are small for S- and P-waves but are significant in the D-waves. 
Agreement with the experimental data is of comparable quality than in the model ig- 
noring vector modes for the low constituent quark mass case. This conclusion has been 
sharpened further by studying the scalar pion radius. Only for low constituent masses 
one can describe this quantity within the Hartree approximation of the (extended) NJL 



model pi. Such low constituent quark masses are also obtained in the currently popu- 



lar estimates of the chiral perturbation theory low-energy constants from extended NJL 
models Further research in this direction should be concerned with a consistent 



implementation of pion loop effects (for some first attempts see e.g. ]3(J ) 
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APPENDIX 



Here we give expressions for the main low-energy characteristics of the tttt scattering 
derived in the 0(q A ) approximation from the amplitude (105) + (111). They coincide 
with the results of |20[ . 
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The quantity Z3 appears in the chiral expansion for the pion mass 
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For the pion decay constant one can obtain 
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Figure captions 

Fig.l. One-loop quark diagrams corresponding to the amplitude A(s, t, u). (a) p meson 
exchange; (b) scalar particle exchange; (c) other contributions. The structure of meson 
vertices is shown in other figures. 

Fig.2. Vertices p — > tctt and o — > tttt in the model with linear realization of chiral 
symmetry. 

Fig.3. Diagrams corresponding to vertex (c) (see Fig.l) in the model with linear real- 
ization of chiral symmetry. 

Fig. 4. Vertices p — ► Tin and a — > tttt in the model with nonlinear realization of chiral 
symmetry. Added to these diagrams should be effects of na\ mixing on lines with 

Fig. 5. Diagrams corresponding to vertex (c) (see Fig.l) in the model with nonlinear 
realization of chiral symmetry. Added to these diagrams should be effects of irai mixing 
on lines with vertex £ M (with the one pion field only). 

Fig. 6a. The S-wave phase shifts for the heavy mass case (in degrees). Upper panel: Sq, 
lower panel: Sq. Data for Sq are from Refs.^TJ and for Sq from Refs. fl&ql . 

Fi g. 6b . The P-wave phase shift for the heavy mass case (in degrees). Data are from 
Refs.|3l. 



Fig. 6c. The S-wave phase shifts for the low mass case (in degress). Upper panel: Sq, 
lower panel: Sq. 

Fig. 7a. The full scalar quark form factor, T, is the sum of the bare coupling of the 
operator m(uu + dd) (double line) to the quarks and a contribution from rescattering of 
the quarks through the four-fermion interaction in the scalar channel, with strength Gs- 
The latter can be represented as a coupling to the composite meson a. 

Fig. 7b. The scalar pion form factor obtained by attaching two pion legs to V. 
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Figure 3 
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Figure 4 
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Figure 6b 
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Figure 6c 
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Figure 7 



39 



